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Abstract: We study the effects of non-trivial initial quantum states for inflationary fluctua- 
tions within the context of the effective field theory for inflation constructed by Cheung et al. 
which allows us to discriminate between different initial states in a model-independent way. 
We develop a Green's function/path integral based formulation that incorporates initial state 
effects and use it to address questions such as how state-dependent is the consistency relation 
for the bispectrum, how many e-folds beyond the minimum required to solve the cosmological 
fine tunings of the big bang are we allowed so that some information from the initial state 
survives until late times, among others. We find that the so-called consistency condition 
relating the local limit of the bispectrum and the slow-roll parameter is a state-dependent 
statement that can be avoided for physically consistent initial states either with or without 
initial non-Gaussianities. 
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1 Introduction 

Both theoretical [1-3] as well as observational studies [4, 5] of non-Gaussianity in the statistics 
of metric perturbations have the potential to open up a veritable treasure trove of insights 
concerning the nature and evolution of quantum fluctuations generated during inflation. 

One of the basic questions about quantum fluctuations of the inflaton concerns their 
initial state. The usual logic for determining this state makes use of the fact that at short 
enough distances or equivalently, early enough (conformal) times, the space-time appears 
approximately flat. Then we can choose the linear combination of solutions to the full mode 
equation that matches to positive frequency flat space modes. This choice leads to the usual 
Bunch-Davies vacuum [6]. While this procedure seems reasonable on the face of it, its basic 
premise is that the inflaton field will exist as a fundamental degree of freedom down to 
arbitrarily short distances/early times. This is a radical assumption; it is much more likely 
that the description of inflation as being driven by a scalar field is an effective one, valid only 
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up to an energy scale A^,. If this is indeed the case, the choice of initial state becomes a far 
more complex issue [7-9]. 

There have been a number of works [10-16] attempting to use the so-called shape of 
the bispectrum, i.e. whether the three-point function of the perturbation is peaked in 
momentum space when the triangle the three momenta form is squeezed, equilateral, flattened, 
or orthogonal (peaked on both equilateral and flattened triangles), to constrain models of 
inflation. More recently, the effects of more general initial states of inflaton perturbations on 
cosmological observables has been studied [12, 13, 17-21], with a particular emphasis on the 
halo bias [22-24]. 

Our aim in the present work is to make use of the recently developed effective field theory 
(EFT) of inflation [25], which incorporates the interpretation of metric perturbations as the 
near Goldstone mode of spontaneously broken time translations together with an effective 
description of the initial state for this mode. Doing this naturally incorporates the fact that 
there may be a limit to the domain of validity of this description via a cutoff indicating 
when the Goldstone theory becomes strongly coupled. We develop a method for integrating 
the initial density matrix into the effective action within the context of the in-in formalism 
[26-29]; this allows us to use methods of in-in perturbation theory in evaluating the effects 
of non-trivial initial states on correlation functions of the metric perturbations. 

This combination of formalisms allows us to answer a variety of questions, such as whether 
the so-called consistency relation [30, 31] for the bispectrum in the squeezed limit is true for 
an arbitrary state. We find that within the context of the EFT of inflation, non Bunch-Davies 
initial states can in fact violate this condition allowing for a loophole in the claim that the 
observation of local non-Gaussianity would necessarily rule out single field inflation. 

In section 2 we review the in-in formalism for general initial states. We take particular 
care in showing how to construct Green's functions that incorporate the choice of a non- 
standard initial state. Section 3 then reviews the EFT of inflation described in [25] . Section 
4 uses these techniques to impose constraints that the initial state must satisfy in order 
to be a state consistent with an inflationary phase as well as with the painstaking WMAP 
observations of the power spectrum. 

Finally, in section 5 we turn to the calculation of the bispectrum in the presence of 
both, the cubic operators appearing in the EFT action for the Goldstone mode describing 
the fluctuations, as well as those from initial state non-Gaussianities. We then conclude in 
section 6 with some further directions this work could be used for. 

2 The in-in formalism for general initial states 

In order to evaluate time-dependent expectation values of operators as opposed to S-matrix 
elements connecting in and out states, we use the in-in or closed time path formalism [26-29]. 
We review this formalism here, with an emphasis on the dependence of these correlators on 
the initial state. We perform most of the analysis here for a scalar field in flat space for 
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simplicity; the conformal flatness of Friedmann-Robertson- Walker (FRW) spacetimes makes 
the translation from the flat spacetime case to the FRW one rather simple. 

Suppose we want to evaluate the expectation value of an operator 0{t) at a time t in the 
Schrodinger picture. This is given by 

{0){t)^TT{p{t)0{t)), (2.1) 

where p{t) denotes the density matrix of the system involved, satisfying the Liouville equation 

.dp{t) 



dt 



[H,p{t)], (2.2) 



with H the full Hamiltonian of the system. If we define the time evolution operator U{t,tQ) 
satisfying the Schrodinger equation 

.mt^ = HU{t, to), U{to,to) = X, (2.3) 

we can solve eq.(2.2) as 

p{t) = U{t,to)pito)U\t,to). (2.4) 
Inserting this into eq.(2.1) we then arrive at 

{0){t) = Tr {uHt,to) 0{t) U{tM p(io)) , (2.5) 

which can be read as follows: start with the initial state pito) and then evolve it to a time t 
at which point the operator 0{t) is inserted. Finally, evolve back to ^o- This is the origin of 
the closed time path required to describe in-in expectation values. 

We can arrive at a more useful formulation of this result by explicitly evaluating the trace 
in the field representation where we use the eigenstates of the field operator ^>(x, t), 

^{x,t)\M-)) =Umt{-))- (2.6) 

We can write the expectation value now as 

Tr(c/t(t,to) 0{t) U{t,to) p{to) 

d^t, d4 d(j)t d4>t {^to\uKt,to)\cl)t) {cl)t\0{t)\4>t) {4>t\U{t,to)\4^t,) {4>to\p{to)\<Pto)- 

(2.7) 

Using the standard representation of matrix elements of the time evolution operator in terms 
of the path integral we can write the above expression as 

{0)it)= [ P<D+ P<D-exp(i (5 [$+] -S [$-])) 0(t) p (<!>+, <I>-;io), (2.8) 

Jbc 
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where the fields in p (<&^, <I>~; to) are evaluated at the initial time to and 5 [$^] is the action 
of the system. The subscript be indicates that the path integrals are taken only over the 
configurations satisfying ^~^{x,t) = <l>~(x, t), ensuring continuity of the field viewed as being 
defined over the entire closed contour. 

We can abstract a generating functional Z [J^, J~;to] from eq.(2.8), 

Z[j+,J-;to] = yp1>+ Pci>-exp(i(5[cD+,J+] -5[c|>-,J-])) p («3>+, to) , (2.9) 

where we have added sources to the appropriate actions. We can also extend the contour 
to t — 7- +00 by turning the sources off after the latest time in the string of field operators we 
may be computing the correlators of. 

This generating functional can be used to construct a perturbative expansion that takes 
the existence of non-trivial initial conditions into account. In the Appendix, we show how 
to derive the in-in Green's functions for an arbitrary Gaussian initial state; initial state non- 
Gaussianity can then be treated as an interaction term with support only at the initial time. 

For use in the next few sections, we give here the parameterization of the initial density 
matrix, 

p{^+,^~-M) = iVexp(icS[<^+,<I>-;to]), (2.10) 

with 

S [^'+,$-;to] = 52 to] + 5int to] , (2.11) 

N being the normalization chosen so that Tr [p[tQ)) = 1 and the subscript int indicating higher 
order interactions in the initial state. The action 5 is in general complex, with hermiticity 
of the density matrix imposing the condition 5 [$"*", to]* = —5 [<I>^, <I>^; to]. We consider 
the quadratic and cubic part of what can be thought of as a boundary action on the initial 
time hypersurface, 

52 [^>+,$-;to] = ^ y d^xd^y {^+{xM)A{x-y)^+{y,to)-^-{x,to)A*{x-y)^~{y,tQ) 

+ 2i^+{xM)B{x-y)^-{y,to)], (2.12) 
53[$+,$-;to] = W^^xi {C(xi,f2,:?3)$+(xi,to)$+(f2,to)^+(x3,to) 

i 

- C*(xi, X2, X3)$~(xi, to)«>~(x2, to)^~(x3, to)} 

^\ j W^^^i {L'(xi,f2,X3)$+(xi,to)$+(x2,to)$"(f3,to) 

i 

- i?*(xi,X2,X3)$~(xi,to)^'-(x2,to)«>+(x3,to)}, (2.13) 

where we used spatial homogeneity explicitly in the quadratic part; it is implied in the cubic 
action. The kernel C{xi,X2-,x^) is fully symmetric in its spatial arguments, while D{xi,X2,X3) 
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is only necessarily symmetric in its first two. Note that all kernels are defined at the initial 
time tQ. Further the kernels can include non-local interactions, which is not in contradiction 
with the fact that the field theory is local. For future reference, we rewrite the above terms 
in momentum space. 



to 



1 

3! 



n 



+ i$i(to)Sfc$:.(to) + i<l>z{to)Bk<^+^{to)], (2.14) 
d^A;. 



(2.15) 



(2vr)3 

X {%,.^..^3'^J(*o)'l>t (to)<^t (to) 

with a similar term corresponding to the kernel D^^ in the cubic action. 
3 The effective field theory of infiation 

The EFT developed by Cheung et al. [25] mimics the construction in spontaneously broken 
gauge theories. At sufficiently high energies, the dynamics of the longitudinal gauge degree of 
freedom can be described by that of the would-be Goldstone boson. In the cosmological case, 
instead of an internal gauge symmetry, time reparameterization symmetry is spontaneously 
broken by the choice of inflaton zero mode configurations (/'o(i); inflaton fluctuations TT{t,x) 
are then defined via 



<P{t, x) = (^o(i) + Hit, x) = Mt) + (Mt + vr(t, x)) - Mt)) 



(3.1) 



Time reparameterization invariance is realized non-linearly on the ir field; if t — t- = t+^^{t, x) 
is an infinitesimal time diffeomorphism, then the linear realization on (j){t, x) is now enforced 
non-linearly on 7r(t, x) through 



t + TT{t,x) = t' + TT'{t',x) 



(3.2) 



The logic is to first choose a gauge for which vr = and construct the most general form 
of the action, which is invariant under time- dependent spatial diffeomorphisms organized as 
fluctuations about a background solution. 



S 



I 



d^x 



1 nji 
-g -Af2^i?-A(t)-c(t)/'' + 5]-ji(5™ + l)" + <5ir,, terms , (3.3) 

n=2 

where we use the FRW metric, g^i, = (— 1, a^(t), a^(t), a^(t)). For simplicity, we neglect 
functions of the perturbed extrinsic curvature since, in the simplest scenarios, these are higher 
order (see [25] for a discussion of exceptions to this) . This form is clearly sufficient to capture 
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the p{X,(f)) models, X = —-^g^^d^cfiducf) being the canonical kinetic term, for instance. By 
re-introducing the time reparameterization gauge transformation a la Stiickelberg, the action 
can be made fully diffeomorphism invariant, albeit in a non-manifest way, as follows. 



S 



d^x y/—g 



Mi 



Ml^R - A{t + it)- c{t + it) ((1 + + 2(1 + ir)diTTg'^' + g'^diTrdjTr) 



+ X] ^ (1 + (1 + ^) V° + 2(1 + ir)d,7Tg'>' + g'^ d^Trdjir)'' + 



n=2 



(3.4) 



with overdots denoting cosmic time derivatives. The functions A(t) and c(t) are determined 
by the requirement that the one-point function for vr vanishes. At the classical tree level this 
implies 



H 



1 



3M2j 



Ml, 



[Ctree(t) + Atrce(t)] , 
■ Ctrcc 



(3.5) 
(3.6) 



H = a/a being the Hubble parameter. However, at the quantum level, for a given quantum 
state, there will be tadpole corrections. They encode how a given quantum state backreacts 
onto its background. For now we will keep A(t) and c(t) arbitrary, anticipating that both will 
receive quantum corrections. 

The clear advantage of the EFT approach comes from recognizing that physics at higher 
energies is dominated by the dynamics of vr. In the high energy limit we can neglect fluctua- 
tions in the geometry, and neglect all non-derivative suppressed interactions of vr (except the 
linear ones which are necessary to determine the tadpole cancellation condition), giving in 
this limit the following effective action. 



Sn = j d^x a^{t) 



a{t)'K -\- c{t) I 7T 



•2 



a2(t) 



n=2 



^2 



aHt) 



(diTT) 



where we have integrated by parts and defined 

a{t) = A + C + 6Hc. 



(3.7) 



(3.8) 



This action gives a general EFT for inflation, from which we can recover specific models by a 
judicious choice of parameters. At cubic and even quadratic order it is also possible to have 
additional terms with higher spatial derivatives, as considered, for example, in [15]. However, 
these terms are always expected to be subdominant to the above ones, unless the coefficients 
of the terms in eq.(3.7) happen to be unnaturally small. 
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The restriction to high energies means that the decouphng effective theory will break 
down shortly after horizon crossing. However, the simple relation tt = —Q/H between the 
Goldstone mode vr and the curvature perturbation ^, valid in single field inflation, tells us 
that we can use the EFT action to follow a given vr mode up to and just beyond horizon 
crossing, and then use the constancy of C, outside the horizon [25]. 

In order to compute correlation functions such as the power spectrum and the bispectrum, 
we can truncate the action, keeping only terms out to cubic order in vr and its derivatives. 
We also remove the tadpole term so as to ensure that quantum corrections to the slow-roll 
parameters are small and the action is at least quadratic in vr, as discussed in the next section. 
Moreover we include leading contributions from the mixing of vr with gravity, by expanding 
the lapse and shift to first order in vr. We can then rewrite the Lagrangian density to cubic 
order following [31] as 

C = Ci+C2 + C^, (3.9) 
where Ci is the tadpole term above while 

C2 = a\t)M^ {t:^ - [dii^f + 2,eH^T^^^ , (3.10) 

£-3 = a'{t) (c^s^'' + ^^^^ i^^7T)^ + C,^2 7r^2 ^ ^^l|^vr {diTrf + Cnl^^.ttS^^tt^ , 

(3.11) 

at next to leading order in slow-roll. The final term is a non-local interaction, which is 
just an artifact of the choice of gauge; all gauge-invariant observables will be local. Here 
= eMpji?^/Cg, e = —H/H^ being the usual slow-roll parameter, and we have defined the 
effective sound speed for perturbations as 



The coefficients of the cubic terms are given by 

C,.=M\l-cl)(^l + ^^y (3.13) 

C,^2 = M^H (-6e + 5-2S + Sec'i - 2e^ (l - c^) j , (3.15) 

C^(9,)2 =M^i/(e-5c2), (3.16) 

Cnl = M^H (^^^ , (3.17) 

where 5 = —e/eH and s = Cs/cgH are other slow-roll parameters. 
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For use in later sections we convert the action to conformal time t] (related to cosmic 
time t via dry = dt/a{t)), with primes denoting conformal time derivatives, 



£2 = a\r])M^ [tt'^ - {d^nf + 3ea\r])H\^j (3.18) 
^3 = a(77) (c,37r'3 + C7^(^^)2vr' (a^vrf 

+ a(??)C^^2 7r7r'^ + a{r])C^^Q^^2TT {di'Kf + a(r/)CNLvr'5i7ra'^7r'^ . (3.19) 

The overall factor of a(ry) comes from the change from cosmic to conformal time in the 
measure. 

We can rewrite the Lagrangians above in terms of a field x with the correct mass dimen- 
sion and kinetic term. Thus define 

X{x,v) = V2M'^a{r])-K{x,r]). (3.20) 

Then we can write, with Ji = a' /a, 

CI = \ (x'^ - cl [d^xf + + 3e^^) x^ , (3.21) 
\a{7])M^ a{r])M ) ^ y a{7])M^ a{r])M j ^^'''^> 

+ ^ {x - nx) dad'^ [x! - nx) , (3.22) 

at next to leading order in slow-roll. We have defined 

C^3 _ C^3 ^ir{dnf _ ^nidnf 



2V2M6 A?2' 2V2M6 M2 ' 

2V2Af6 " ~' 2^/2M6 = M ' 2VM^ ^ 



(3.23) 



where the coefficients with tildes are dimensionless. 

This is an effective field theory and as such, has a limited domain of validity. There are 
two important scales for the EFT [16]; the symmetry breaking scale at which time trans- 
lations are spontaneously broken by the background evolution and a description in terms of a 
Goldstone boson first becomes applicable and the strong coupling scale A^, where perturbative 
unitarity is lost. The strong coupling scale can be computed [16, 25] and is given by 



A^ = 2eM^yH^r^ = i ] M\ (3.24) 



,5 



2cl 



-8- 



We can use the existence of the scale A^, to define what we mean by the initial (conformal) 
time t]q which we will use in computing the relevant correlation functions. We define it to 
be the earliest time at which the EFT is valid for scales that appear on the CMB sky today. 
While inflation could have started much earlier, a description of the physics of the scales we 
measure today would require strong coupling physics at those earlier times. 

4 Constraints on general initial states 

Our formalism allows for the use of any state as an initial state. However, not all states are 
physically acceptable. There are two types of constraints that must be imposed on the vr, or 
equivalently the x theory. The first type is a consistency condition due to the fact that the vr 
theory is an effective one and all of the physics we use this theory for must remain consistent 
with the precepts of effective field theory. In particular, we should not be able to excite modes 
near the cutoff A*. One way to do this is to absolutely forbid that such modes appear in the 
state, and this is the approach we will take here. It's worth noting, though, that this may 
be too strict. We could imagine a scenario where we have some particles with energies near 
the cutoff in the initial state, but their contribution to the energy density is small and the 
time evolution of the system is such that we do not produce any more such particles. The 
exponential expansion associated with inflation only helps to enforce this weaker criterion. 
This may be worth pursuing in future work. 

For now, we will assume that the state satisfies the Hadamard condition [32], which 
demands that the Bogoliubov /J^ coefficients fall off faster than k~'^ at large k. As in [13], 
we enforce this by demanding that — )• for /c > A^ao/cg (the factor of gq = ^(r/o) comes 
from choosing A^, as the physical cutoff), that is to say, as the mode momentum is redshifted 
across the cutoff of the theory at i]o, it starts off in its vacuum state. 

Beyond this, the next set of constraints enforces the fact that inflation should have 
occurred. This takes two forms; the backreaction of the energy density in the initial state 
should be subdominant relative to Mp^H^ and corrections to the slow-roll equations of motion 
should be much smaller than the scale of the original terms themselves. We will take these 
up in turn. 

4.1 Backreaction 

To impose the backreaction constraint we demand that p{r]) <^ Mp^H'^, where the energy 
density p is given by {T^q}- This will generally involve divergent mode integrals and a more 
sophisticated analysis would involve using something like an adiabatic expansion [33] of the 
modes to isolate the divergent terms. The higher dimension terms in the action would then 
induce divergences beyond those of the free theory, but these would be absorbed in higher 
dimensional counter terms in the usual way when dealing with effective field theories. To arrive 
at the estimates we are looking for, we will default to a cruder method where we compute the 
expectation value of various operators, subtract off their values in the Bunch-Davies vacuum 
and then cut the integrals off at Cgk = aoA*. 
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The quadratic part of the Lagrangian will give rise to contributions to the energy 
density such as 

When the X^iv) expanded in terms of the Bunch-Davies modes hj^{rj), there will be terms 

proportional to , ^^(^) > well as cross-terms proportional to h'i{r]), h'?{ri) with 

their complex conjugates. We want to argue that we can neglect these latter terms. They 
are proportional to exp {±2icskr]) and we would like to be able to use the Riemann-Lebesgue 
lemma to show that the oscillatory nature of these terms will wash out the integrals. This 
requires a large parameter in the exponential; at early times we can replace r] with 770 and 
ferret out this large parameter by writing 

^ aoH oqA* H 

Once the cutoff A^, becomes smaller than H, the use of the EFT to describe inflation is no 
longer valid, so we expect that K^/H ^ 1. At late times, on the other hand, the factor of 
a~^{r]) in the energy density will exponentially redshift the cross-terms. 

As an example consider the contribution from the spatial gradient part of the kinetic 
term, 



1 f d^k 



aHv)J (27r 



\3 



elk' (2\p>\' + l) \h^M'^ (4-3) 



where we use the notation in the Appendix for the modes. We should note that the above 
result is for a pure state; for a mixed initial state, we need only multiply the integrand by 
(Tfc (see eq.(A.45)). We can replace ^2 + 1^—7-2 |/3^|^, likewise for a mixed state we 

can replace cjfe ^2 |/3^ |^ + 1^ — )• 2ak + (cfc — 1)) since the remainder is the Bunch-Davies 
contribution we are subtracting, and write the modes in the de Sitter limit, 

(1.4) 

We also model as (3^ = f{q), where q = Cgk/ (aoA*) and f{q) falls to zero faster than 
at high q. Putting all of this together we find a constraint of the form 
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with K = A^/H ^ 1 and q„im ^ 1 an infrared cutoff. The left hand side can be bounded 
above by 

tj2 \2 roo 

Hf/ d,,A,)(i + .V) 

lj2 \2 ( rl poo ^ 

= / d,+ / dg (,A,)(1 + kV)). (4.6) 

For the second integral we can use the fact that f{q) ~ fo/q^'^^ , > 0, while we can model 
the result of the first integral by saying that f{q) = /o, (?min < ^ < 1 to arrive at the result 
that the integral can be approximated as J'^k'^ up to factors of order unity. We can then 
satisfy the backreaction constraint by demanding that 



But the left hand side can be rewritten as /gC^Mpjif^e/ (vr^ (l — c^)) so that the constraint 
becomes f'^ <C vr^ (c^^ _ l) /e. Of course, this will change if different models for the Bogoli- 
ubov coefficient are chosen, but we can see that the constraint is relatively easy to satisfy 
even with /o of order unity. 

For a mixed initial state, we can model {at — 1) as (crfc — 1) = s(g), with s{q) decaying 
faster than q~^ at high q. Then taking s{q) = sq (with sq positive) in ^min < Q < 1, the 
above constraint becomes /g (1 + sq) + sq/2 ^ vr^ (c~'^ — l) /e. We see that this constraint 
can again be satisfied with both /g and sq of order unity. 

The cubic interaction terms in £3 also contribute to the energy density. While we would 
have to go to higher order in cubic interactions to find an effect within the Gaussian part of 
the initial density matrix, we can find a first order contribution from these terms when we 
include initial non-Gaussianity, as encoded in the cubic action of eq.(2.15). In fact we can 
use cubic interactions to place bounds on the high behavior of the kernel Cg^ g^. Let's 
show how this works by considering the contribution to the energy density of the operator 
x'^ appearing in the cubic Lagrangian of eq.(3.22). We have 

^Pcubic(^) = §^^(x'^(x,r?)> (4.8) 

= I n0^"^''^'^ {^,M^,MX,M ■ (4-9) 

We then use the cubic part of the initial state action to compute the momentum space 
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expectation value as 



^ ( ^kiMM (^)^k ^''o^^J ^''o^^J ^'^o) 



Gaussian 



(4.10) 



The expectation value can be computed in the Gaussian theory by taking the time i] to lie on 
the + part of the time contour and then using the Green's functions found in the Appendix. 
Using the symmetry of the kernel Cg^ in its arguments we find for example that 



Xffi {'n)x^'{'n)x]f'{ri)xl {m)x^ (^o)xi' (^o) , ^ . 

^2 m ki fe2 ks I Gaussian 



% 

Z\'W{2'KfbHk^i) a/,>(r?)/,<(%)). (4.11) 



Inserting this back into eq.(4.9) and doing the qi integrals we see that the cubic contribution 
to the energy density is 



1 



M2 a^{ri) 



n 



{27TY 



(4.12) 



The strongest backreaction constraint on the kernel C^^ will come when we only keep 
the terms in the mode functions containing the Bogoliubov coefficients, since these are 
order unity. Doing this we can bound the cubic contribution as 



5pcubic(??) ^2%^ 



n 



< 



< 



M2 ae(rj) J ^ (2^)^ 



(2vr)3 6' (Y^ k 



C. 



ki,k2,ks 



{drihk,{v)h*k,ivo)) 



1 



ao 



167r4c6 \a{7]) ) M'^ ^ 
1 A! 



dqi d(?2 ql ql 



C, 



levr^cf M2 



C 



ki,k2,-{ki+k2) 



(4.13) 



where we have used the de Sitter form of the Bunch-Davies modes hk{if) and kept only 
the leading terms in the integral. We have also assumed that Cg^ .^^g^+fc^) off faster 
than at high k so that the integral from g = 1 to oo is subdominant, and treated 

^fci k2 -{ki+k2) essentially constant in ^niin ^ ^ ^ !• The actual bound one would get will 
of course depend on this angular dependence, but we just want to get a rough estimate as to 
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how big the kernel can be. Enforcing that this last result be much less than Mp^H^, using 
the expression for A^, and taking C^3 ~ 1 yields the bound 



« ^ i'-'sY ■ (4-14) 



In the small Cg limit, this is not a very stringent bound. A similar bound on Dj^^ obtains. 

The point here is not the exact bound (we have been somewhat sloppy in our estimates); 
rather it is the fact that we can get some constraints on the initial non-Gaussianity from 
backreaction. In the next subsection we will consider corrections to the background equations 
of motion and get further bounds from that. 

4.2 Tadpole 

The next constraint has to do with not only ensuring that the energy density associated with 
a non-vacuum initial state does not disrupt the inflationary era, but that the inflaton field 
zero mode actually follows the appropriate equation of motion. This is equivalent to the 
requirement that the tadpole of fluctuations about the correct zero mode trajectory vanishes. 
Going back to eqs.(3.7, 3.8), we see that the coefficient of the terms proportional to TT{x,r]) 
in the action are of order eMp^H^; in terms of the canonically normalized field x these terms 
are of order eMp^H^ / M"^ ~ y/eCsMpiH'^ . Thus, corrections to the tadpole (^X^iv)'^ due to 
changes in the initial state must be much less than this. 

As in the previous subsection, we will compute these corrections both due to the change 
in the modes in the Gaussian part of the theory as well as those due to the introduction of 
non-Gaussianity in the initial state. In computing the former corrections we will follow the 
prescription we used in computing the backreaction, i.e. we will subtract out the Bunch- 
Davies part and then use the fall off of the (3^ Bogoliubov coefficient to cut the momentum 
integral at k = aoA^/c^. The strongest new constraints on will come from the use of 
the dimension 6 operators (when written in terms of x) 

(diir) since their coefficients 

are not slow-roll suppressed. For ease of calculation, we focus on the second of these two 
operators, though the constraint from the first is essentially the same. 

The zeroth order contribution to (x'^iv)) coming from the linear term in eq.(3.7) is 



(xf(^))^°^ = ^{2nf5' {o) j\' a\r,') (G++(r?,r?') - G+-(r?,V)) ("^^j • (4-15) 

Written in this form, the correction to the tadpole due to the new modes can be written as 
the following correction to 0(77), 

= ^2c,„,, / (ai \JSm + n\ISM-) ("6) 



If we use the de Sitter form of the modes, we find for a pure initial state 




+ 7^^^ = (2 \P>\' + 1) ( cM - \\ ■ (4.17) 
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Subtracting the Bunch-Davies contribution, using the high-A; fall-off properties of , and 
again transforming to the variable q = Cgk/ (qqA^) allows us to write 



Here we have modeled the Bogoliubov coefficient in the same manner as for the backreaction 
calculation of the previous subsection. The tadpole constraint is then satisfied by taking 
/o <C 2\/27r2 ~ 10 for small Cg and C'^^g^^2 ~ 1. Note that this constraint is parametrically 
smaller in e than the one coming from backreaction; however, it is not itself parametrically 
small in e. As earlier, we can generalize this constraint to the mixed initial state case: 
/2 (1 + so) + so/2 « 2V27r2 ~ 10. 

Next we calculate the contribution of non-Gaussianity in the initial state to the tadpole. 
We find 



6a{v) = -M' I ^6{^ - ,?o) / jz-sC^,,.n,6 ^^^^^ 1 ' (4-19) 



V2ao J (27r) 







where 6{rj — tjq) comes from writing G~^^{rj, rjo) as an integral over r]' so as to match the form 
of the tadpole in eq.(4.15). Once again using the known high-fc behavior of both as well 
°f ^kiM,-{ki+k2) find that 

- 5a(r?) = M \ /I S{7j - rjo) f ^C, (l + , (4.20) 

and hence the tadpole constraint can be written as 

« ATT^cr (1 - c'gf H, (4.21) 

«o 

where we have used rotational invariance to write = Cq and assumed, as in the back- 

reaction case, that it remains essentially constant in the region of integration. To understand 
what the delta function implies, integrate both sides with respect to r] from 770 to r/ = to 
see that 

Cq^i « 8^2^3/2(1 _ ,2)1/2 ^_a,r^^H) = Svr^c^/^ (l - c',f' . (4.22) 

Note that this is a constraint on the (extremely) squeezed limit of initial non-Gaussianity and 
that it, just like its Gaussian counterpart, is parametrically smaller in e than the backreaction 
bound, though it itself is not e suppressed. 
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4.3 Power spectrum 



Our current knowledge of the power spectrum of fluctuations, in particular the fact that it 
is nearly scale-invariant, can also serve to constrain the form of the initial state. The power 
spectrum Pk for the curvature perturbation (" is defined in terms of the two-point function. 



fci ^te 



The dimensionless power spectrum Vk = ^^-Pfc is usually parameterized in terms of the 
measured normalization A{kp) and spectral tilt Hs as 

(4.24) 



(4.23) 



V^ik) = A{k,) ( - 



kj^ 

kp being the pivot scale fixed at 0.002 Mpc~^. The 7- year WMAP observations combined 
with BAO and Hq measurements [4] fix these parameters at 



A{kp) = (2.430 ± 0.091) x 10^^ 



0.968 ±0.012, 



(4.25) 



at the 68% confidence level. The power spectrum generated from general initial states for the 
perturbations must be consistent with these constraints. 

Using eq.(A.48) and the expressions for the mode functions, we find that the power 
spectrum for general initial states is given by 



Vdk) 



k^ 



1 



27r2 2M^^ea^{r]) 
k^ cl 1 



27r2 2M2 ea2(r?) 



87r2M2jeCs 



(Tk \al 



Pk 



(4.26) 



Csk=aH 



where again is the mixing parameter defined in the Appendix and a^, are the Bo- 
goliubov coefficients. We are allowed to calculate observables at late times [i] — )■ 0~) since 
the perturbations freeze out at horizon crossing {cgk = aH). Here we have only calculated 
the leading contribution to the two-point function since, in the absence of an initial state 
trispectrum, the next contribution is an eight-point function, with three additional x fields 
coming from initial state non-Gaussianity and another three from cubic interactions in the 
EFT action. 

The constraint that the power spectrum is nearly scale-invariant translates into the con- 
dition that 

d 



d In k 



In 



Hi 



< 



0(0.01). 



(4.27) 



Csk=aH 



which is similar to that found in [22]. How much of a constraint this is for the initial state 



will depend on the behavior of the P?' coefficients for k's near horizon crossing. 
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5 Calculating the bispectrum 



While measurements of the power spectrum serve to place important constraints on allowed 
initial conditions, it is basically a measurement of only two numbers: the amplitude and 
the spectral tilt. The bispectrum or three-point function promises a much greater amount 
of information seeing as it is a function of the magnitudes of, and relative angles between, 
three vectors forming a triangle. The bispectrum i?^^ for the curvature perturbation 
is defined as 

The bispectrum directly probes the dynamics and interactions of the inflaton, and so we 
expect it to be very sensitive to the initial state of the fluctuations. Different inflationary 
models can be identified by the 'shape' of the resulting bispectrum, defined by the type 
of triangle for which the amplitude of the three-point function is largest. The most well- 
studied shape of the primordial bispectrum is the 'local' type, in which the perturbation (or 
equivalently the primordial gravitational potential) is a simple non-linear function of the local 
value of a Gaussian field. The non-linearity is usually parameterized in terms of a constant 
/nl, 

~ ^fc.Gaussian ^-^^^ (''fe.Gaussian {^fc.Gaussian)) ' (^-2) 

where the subscript indicates a Gaussian random field. The above parameterization generates 
a bispectrum of the form [2, 34, 35] 

^tt,ks = ^^NL {Pk,Pk, + 2 perm.) , (5.3) 

which is maximized in the squeezed limit, in which one of the momenta is much smaller than 
the other two (k-^ ^ A;i ~ k2). Other commonly studied shapes of the bispectrum include 
the equilateral {ki ^2 ~ k^), flattened {k^ ~ /ci + ^2), and orthogonal (peaked on both 
equilateral and flattened triangles) configurations. 

General initial states have been found to introduce interesting features in the primordial 
bispectrum [12, 13, 17-20], including enhancements in the squeezed and flattened momentum 
configurations. The bispectrum for the curvature perturbation can be obtained from the 
three-point function for the tt field, 

(5.4) 

with 53 given by the cubic terms appearing in the version of eq.(2.15) for the canonically 
normalized field x- Our goal in this section is to calculate the effective /nl for local config- 
urations using generalized initial states, with an eye towards understanding to what extent 



- 16 - 



initial state effects can alter the inferences we might make about inflationary physics if and 
when non-Gaussianity is detected. In particular, we ask whether or not the statement of 
the so-called consistency condition, that is, that a local non-Gaussianity in single field infla- 
tion should be parametrically slow-roll suppressed and hence undetectable, still holds within 
the context of the EFT approach to both the interactions as well as the initial state of the 
fluctuations. 

5.1 General Gaussian initial states 

We start by considering the bispectrum for a general Gaussian initial state at leading and 
next to leading order in slow-roll. For Bunch-Davies initial states, it was shown in [31] that 
leading order terms are subdominant in the squeezed limit while next to leading order terms 
are required in order to recover the usual consistency relation of [30]. 

5.1.1 Bispectrum at leading order in slow-roll 

At leading order, = —Hit and the three-point function for ^ is just —H^ times the three- 
point function for vr. For the two leading order operators vr'^ and 7r'((9j7r)^ in the cubic EFT 
action we find that 



A:. (1 



X / dr]' a{rj') 



vo 



1 + 



2 Ml 



a,,G^;++(o,V)) {d,,Gi;^{o,v')) {d,,Gi;^{o,v': 



3M| 



+ permutations -|- c.c. 



(5.5) 
(5.6) 



where a(r/') = —l/{Hr]') and in the last line we have defined the contribution to the three- 
point function from the operators vr'^ and 7r'(5j7r)^. The dot product A;2-^3 here can be 
simplified in terms of magnitudes of the momenta using Y^ki = 0, from which it follows that 
^2-^3 = ^ {kf — k2 — k"^). The -|— |- Green's function for the vr field at conformal times rj > rjo 
is given by (see the Appendix) 



(0,r?') 



1 



4M^,_ff2e a{r])a{ri' 



afc(l - icskr,')e''^'''>' + bk{l + icskr]')e 



-icskr)' 



(5.7) 



with 



AM^^ek 



(5.8) 
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where we have introduced the functions 



1 



- [{a, + 1) (a> - /3>) a>* - {a, - 1) (a>* - (3>*) /3>] 



- [- (a, + 1) (a> - /3>) + i^k - 1) " «fc ] ' 



(5.9) 
(5.10) 



with a^, being the Bogohubov coefficients. Strictly speaking, the derivative of the Green's 
function in eq.(5.8) should have another term coming from the 0(r/o — rj) term in eq.(A.36). 
However, we can imagine calculating the integral in eq.(5.5) from r/ = r/o + e to 0, at which 



point we can ignore this extra contribution. Using the fact that 
rewrite a/;, 6fc as 



ak 
h 



1, we can 

(5.11) 
(5.12) 



It's worth noting that the backreaction and tadpole constraints give bounds on the combi- 
nation of + {ck — Furthermore for o"/; > 1 both and have a term that is 

independent of (5^ , though the constraint on this term is of the same order as that on | f3^ | ^ . 

We are now left to evaluate the integrals in eq.(5.5). On calculating these integrals we 
find that the two terms in eq.(5.6) are given by 



|(2vr)3 53 cs{l 



2 Ml 

1 -I 3 

^ 3 M| 



1 



X 



and 



/,m,n=0 
+ C.C. 



C^C^C^) (??) 



(5.13) 



X 



E 



(/) im) in) 



1 



-k2 ^fcs' ^^'{d^) 



M^je2 k^k2k3 



l,m,n=0 
+ C.C, 



(5.14) 



where 





1. 



(5.15) 



Here J^^n and T^i{J)^y denote the shape functions. 



-^7r'3(Pl,P2,P3,%) 



JCsKirio 



+ 



(. 



2ir]o 



4Kf CsKi \clKl CsK 



(5.16) 
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and 





^icsKiT^o r^6 _ sKfK^ + llKfKl - 4KfK^ - AKiK^Ki + 12KI 



CsKiKl 



(5.17) 



where for brevity of notation we have suppressed the exphcit momentum dependence of the 
functions Ki, K2, and K3, 



Equations (5.13) and (5.14) give the three-point function at leading order in slow-roll. The 
difference from the Bunch-Davies limit is the presence of both positive and negative momenta, 
±fcj, the consequences of which we will discuss below. 

Let us first consider the Bunch-Davies limit of the three-point function at leading order 
obtained here. In this limit we can set r/o — — oo(l it ie) (choosing the plus sign for positive 
frequency modes and the minus sign for negative frequency modes, so that we recover the 
standard vacuum of the interacting theory), and drop all of the exponential terms. With 
(Tfc = 1 and = 1, = we have = 1, 6^ = 0, so that the only term that survives 
in the three-point function is the dfe^a/jjOfca one. In the squeezed configuration the shape 
functions J-^ia and J^^/(g^)2 for this term are proportional to 1/kf, and the full three-point 
function is proportional to 1/ [kfk^^ = 1/ {kfk^) x (k^/ki)"^. Therefore, using eq.(5.3) and 
Pk oc l/k^, we conclude that the leading-order contribution to /nl is suppressed by a factor 
of {k^/ki)^ . This is in agreement with earlier results in the Bunch-Davies limit [31]. 

Away from the Bunch-Davies pure state the story becomes more interesting. The dom- 
inant contribution to the three-point function comes from terms like 0^^6^20/03, which can 
be of order unity. In the squeezed configuration the momentum dependence of the shape 
functions J-^is, and T^i{p^Y for this term is l/Zcg, and that of the full three-point function is 
1/ (A:^A;|) = 1/ (Aif/cg) x (ki/k^). The leading-order contribution to /nl is thus enhanced by 
a factor of ki/k^, showing that general Gaussian initial states can lead to enhancements in 
/nl at leading order in slow-roll. We note that the three-point function is also enhanced for 
the flattened triangle, due to the presence of vanishing denominators in this configuration. 

5.1.2 Bispectrum at next to leading order in slow-roll 

Let us now calculate next to leading order terms in the bispectrum. At this order = 
—Htt + ^eH^n'^ outside the horizon [31]. In Fourier space, = ~^'^k ~^ \^H'^{t^ * 



K2{P1,P2,P3) 
K3{pi,P2,P3) 



Pi +P2 +P3, 
{P1P2+P2P3 +P3P1) 
{PlP2P3f^^- 



1/2 



(5.18) 
(5.19) 
(5.20) 
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where the convolution (vr * vr)^ is given by J ^2^3 '^g-g-^g- next to leading order, all five 
operators, vr'^, 7r'((9j7r)^, vrvr'^, 7r((?i7r)^, and Tr'SivrS^^vr', in the cubic EFT action contribute to 
the bispectrum. We have already calculated the leading order bispectrum from the first two 
operators. Since it is and not vr, that does not evolve outside the horizon, we can calculate 
the contribution of these two operators at next to leading order using = ~jj ~^ ea{rj)H'K 
under the integral. The next to leading order piece of the three-point function is then given 
by 



ki k2 fcs 



(1) 



(5.21) 



where 

ki k2 k'i 
fO 



(1) 



X / dr]' a\r]') 



ci{e-6) Gi;^(0,v') h.ks Gl^^(0,v') Gl'^^M) 



2e 



+ permutations + c.c 



hi k2 I ■/^■ji'' 

+ 



»?->o- 



(5.22) 



(5.23) 



The superscript (1) here denotes that this is the next to leading order contribution to the 
three-point function. 

As before, we can calculate the integrals in eq.(5.22) and obtain the three terms defined 
in eq.(5.23). While the first and third terms are straightforward, the second term is a bit more 
subtle; substituting in for the functions G^'~^~^ {0,r]') and using a(r/') = —l/{Hr]'), this term 
appears to diverge. However, on carefully expanding out the integral, adding the complex 
conjugate for the divergent pieces, and observing that aki = —bki (the subscript / denoting 
the imaginary part), we find that all divergent terms cancel out. The three terms in eq.(5.23) 
are then given by 

1 H"^ 



fci k2 k-i I TT^/: 



32 



(-3e + (5-2s) 



1 



X 



1 

E 4? 4f 4? ((-l)'A;i,(-l)-fe2,(-l)"A:3,r?o 



i,m,n=0 
+ C.C, 



(5.24) 
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(1) 



(e-<5)- 



1 

Z,m, 71=0 



T^^o^y (^{-iykr,{-irk2,{-'^rh,m) 



+ ex., 



(5.25) 



and 



1 



1 



X 



1 



64 ' ' V^ "7 c^M^^ekMz 
(0 » 



E < C 47 -^NL ((-I)'fel,(-l)™fc2,(-l)"fc3,% 
Z,m,n=0 

+ C.C., 



(5.26) 



where, 



2KfKl - KfKl - KiKjKl 
KM 



JCsKirio 



{2KfKl - KfKl - K.KiKi + icsKfKlKlrjo) , (5.27) 



iCsKirjo 



+ 



Kf - 2Kfiq) + Cs {KlKl + kIkI - 2KfK^ - 2KiK^Ki) % 



{KfKl - 2KlKlKl) rii 



(5.28) 



and 



-^nl(pi,P2,P3,?/o) 

piCsKirio 



KfKl 



2Kf - 7KfKl + 17KfKl - AKfK^ + AKiK^Kl 
{2Kf - IKfKl + 17KfK'i - AKfK^ + AKiK^kI) 



-ic. 



Kf {KfKl + KfKl - AKiKl + AK^kI) rio 



(5.29) 



Equations (5.24), (5.25), and (5.26) give the three-point function at next to leading order in 
slow-roll. The difference from the Bunch-Davies limit is again the presence of both positive 
and negative momenta, it/cj. 

In the Bunch-Davies limit = 1 and 6^ = so that the only term that survives is the 
flfciOfc20fc3 one. The squeezed limit three-point function for the operators tttt'^ and 7r(9j7r)^ is 
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proportional to l/{k'lk^), while that for the non-local operator is proportional to 1/ [kfk^j = 
1/ {kfkf) X {k^/ki)'^ and hence suppressed by a factor of {k-^/ki)"^. For general Gaussian 
initial states the dominant contribution comes from terms like akibk2ak3, as in the leading 
order case. The momentum dependence of this term in the three-point function for the 
first two operators is 1/ = 1/ (^1^3) x (^1/^3) and that for the third operator is 

1/ (fcffcl) =1/ (^1^3) X (^3/^1)- So the first two operators give an enhanced /nl signal even 
though the third operator is suppressed. Again, we note that there are enhancements at 
next to leading order in the flattened triangle limit as well, due to the presence of vanishing 
denominators. 

The enhancements at next to leading order in slow-roll we report here match with earlier 
results [12, 13, 18, 20] that used Maldacena's cubic interactions [1]. The enhancements at 
leading order in slow-roll that we find are, however, a unique feature of the EFT approach to 
inflation. 



5.2 Non- Gaussian initial states 



We now add non-Gaussian terms to the initial state. Consider a non-Gaussian initial density 
matrix with 



i 53 



{4} ' {4} ^ ^0] = 2V2M^a3(,o) ^ / n [i^^f (E % 



+ 



(5.30) 



where the pre-factor comes from the conversion between the x and vr fields. This gives rise 
to a contribution 



(nG) 



(5.31) 



to the three-point function, where the superscript (nG) denotes the contribution from initial 
state non-Gaussianities. At late times we find 



(nG) 



X 



C. 



G^;++(o,7?o) Gl';^{o,m) G^;++(o,%) + 



-|- permutations -|- c.c. 



(5.32) 



with G^'^^ {0,riQ) given by eq.(5.7) for r]' = i]q. 

In the absence of an explicit form of the kernel, it is difficult to make any detailed 
statements. However, we can estimate the effect of non-Gaussianities in the initial state on the 
bispectrum today as follows. Suppose we pick an initial state close to the Bunch-Davies one, 
with |r/o| large but finite, and w 1, 6fc 0. Let us further choose the kernel .^^^+^2) 
to be approximately constant in Qyjiin — Q — with q = Cgk/ (tto-^*)i B^nd. vanishingly small cit 
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high k. Then the size of the bispectrum from initial state non-Gaussianities relative to the 
squeezed configuration bispectrum in the Bunch-Davies limit is roughly 



\ BD 



C. 



1 ^/e Afpi ks 



BD 
NL 



H h 



cos {2cskirjQ) . (5.33) 



While this looks somewhat slow-roll suppressed due to the y/e factor, there are a number of 
enhancements associated with this factor. With reasonable numerical values, such as /j 



BD 
NL 



ki,k2,-{ki+k2) 



€ ~ 0.01, Cs ~ 0.1, H ~ 10-^Mpi, and k-s/ki ~ 10"^ this evaluates to 10 
cos(2cs/ci??o)- The kernel is not tightly constrained, as we saw in section 4 so in principle this 
could give a large result. However, since we must average over the cosine of a large argument, 
the effect of initial state non-Gaussianities will be washed out. 

This conclusion can easily be modified if our initial state has ^ for the relevant wave 
numbers; then the cosine factor will involve the small momentum and thus be essentially one. 
This would then lead to a potentially large local non-Gaussianity, obviating the conclusions 
of the consistency relation. 



6 Conclusions and further directions 

The formalism we've developed here creates an amalgam of the EFT for inflationary fluc- 
tuations and the treatment of modiflcations to their initial state as interactions supported 
only at the initial time tjq. Using this we were able to quantify the constraints coming from 
backreaction, the requirement that the correct background held equations were obeyed, and 
the near scale-invariance of the power spectrum. We found that it was possible to satisfy 
these constraints and still obtain significant enhancements to the bispectrum in the squeezed 
and flattened limits. These enhancements appear at leading-order in slow-roll and suggest 
deviations from the consistency relation that are not momentum suppressed or slow-roll sup- 
pressed. 

We were also able to include mixed initial states as well as those supporting non- 
Gaussianities in a relatively simple way that allows for simple calculations of their contri- 
butions to physical observables such as the bispectrum and /nl- 

The effective nature of the tt action requires us to ensure that no quanta of the theory 
are excited to energies greater than the cutoff A* at any time that we require the EFT to 
be valid at. The linearity of quantum theory implies that formally, at least, it should be 
possible to place a quantum system in an arbitrary quantum state that can be built from the 
Hilbert space (i.e. Fock space) of the quantum theory. In the context of EFT, however, this 
is no longer the case. If the EFT has a cutoff A*, states containing particles with energies 
greater than or near A* do not make physical sense. A minor exception to this arises if we 
choose single particle states in a theory in which Lorentz invariance is exact or only partially 
broken, since then we can always boost to a frame for which the energy is below the cutoff. 
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Nevertheless, even in an (almost) Lorentz invariant theory for a typical state which is built 
out of multiple particle states with energies E > A^, there is always some Lorentz invariant 
combination that lies above the cutoff. 

Consider a fluctuation with comoving wave number k corresponding to a physical scale 
on the CMB sky today. This fluctuation mode became amenable to being described by the 
EFT when its physical size became of order A*, i.e. Cgk = OEFT.fc-'^* while it exited the 
inflationary horizon when Cgk = akH. This implies that AA^ = InA^,/// is the number of 
e-folds between when the scale enters the domain of validity of the EFT and when it leaves 
the horizon. Bounding A* above by Mpi and taking H ~ 10~^Mpi would give AA^ < 13; 
typically AA^ will be smaller than this. In particular, we can apply these considerations to 
the scale that enters the domain of validity of the EFT at the initial time rjQ-. CgkM = clqA^, 
where oq = a{r]o). We will want the wavenumbers relevant to the description of the CMB sky 
today to be less than kM so that their physical sizes will be less than A^, at rjQ. This means 
that the weak-coupling description of inflation given by the EFT can be valid for at most 
another AA^ e-folds past the bare minimum required to solve the horizon problem. 

Having constructed this formalism, we can use it to address a number of interesting 
questions. First comes the issue of the consistency condition. For non-Bunch-Davies initial 
states both with or without initial non-Gaussianities, we find that the bispectrum in the 
squeezed limit is enhanced at leading order in slow-roll, leading to violations of the consistency 
condition for single field models of inflation. If this condition is not verified observationally, 
therefore, it may not necessarily imply that single field inflation is ruled out as is often claimed, 
but may only be a signature of a non-trivial initial state for the inflaton. It is important to 
note that the terms that lead to this violation vanish for Cg = 1, and hence it is only single 
field models with non-trivial initial states and Cg ^ 1 that violate consistency. Of course this 
also depends on the extent of departure from the Bunch-Davies vacuum, i.e. on the size of 
the bk coefficients. 

On the other hand there have been a number of works [36-38] arguing that this condi- 
tion, as well as generalizations of it involving higher point correlators, follow from the Ward 
identities of softly broken conformal symmetry. It would be interesting to see how our for- 
malism interacts with this set of ideas. On a more phenomenological level, we would like to 
follow how initial state non-Gaussianities enter into observables such as the halo bias and 
what those measurements might be able to say about the initial state of the inflaton; these 
projects are ongoing and should yield fascinating insights into whatever mechanism sets these 
initial conditions. 
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A Green's functions for general initial states 

As it stands, the generating functional Z[J'^ , J~ ]t{)\ of eq.(2.9) is not yet written in a way 
where we can calculate the correlation functions of interest. This can be dealt with by writing 
the initial density matrix as part of the action 

p(cI>+,$-;to) = A^exp(i5 [$+,«&-; to]). (A.l) 

The normalization N is chosen so that Tr(/3(to)) = 1- The quadratic part of 5 to] 
corresponds to the Gaussian part of the initial state and will be used to incorporate the initial 
state into the Green's functions. Higher order terms will be used to input non-Gaussianity 
into the initial state and can be viewed as interactions. We write 

5[$+,^>-;to] = \j d?x d^y {$+(x,to)^(f-y;to)$+(y,to) 

- $-(f,to)A*(f -y;to)«>~(y,to) 

+ 2i$+(f,to)S(x-y;to)$-(y, to) + •••}, (A.2) 

where • • • refers to cubic and higher terms. We have used translation invariance to write the 
kernels only as functions of the difference x — y; rotational invariance forces the kernels to 
depend only on the magnitude of this difference. Hermiticity of the density matrix further 
requires that B{x — y; to) be real. 

Note that in the absence of the term mixing and the initial density matrix could 
be viewed as supplying a boundary action living on the initial time spatial hypersurface 
[39-41]. However, this interpretation is less tenable for a mixed state. 

Consider a canonical kinetic action in the fields, with some potential y[<^^], 

'1 



^kinctici^"^] - 5'kinctic[4' ] = j d'^X 



^ ( t) - , t) ) - y[$+] + v[<^- 



(A.3) 



We can incorporate the quadratic terms in S <l>~;to] into the kinetic action above, and 
similarly higher order terms in S to] into the full action, by multiplying the kernels 

by various time delta functions, 

5 to] = d^xd^y J dt dt' {<l>+{x,to)A{x - y;t,t')<l>+{y,to) 

- {x,to)A*{x - y;t,t')<i>~ {y,to) 

+ 2i<l>+{x,to)B{x-y;t,t')<^-{y, to) + ■■■}, (A.4) 



with 



Aix - y; t, t') = 2A{x - y; to)6{t - t')6{t - to), (A.5) 
Bix - y; t, t') = 2B{x - y; to)J(t - t')<5(t - to), (A.6) 
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where the factors of 2 multiplying the original kernels are due to the fact that we are taking 
the time integrals to run from t = and that 

dt 6{t - to) = J. (A.7) 

to ^ 

With this parameterization of the initial state in hand, we can now turn to the construction 
of Green's functions that encode the effects of the initial state. 

Using spatial translational invariance, we can consider the kinetic operator as acting on 
the momentum modes (t) of the field. If we define the vector 

the kinetic action including the initial time kernels can be written as 

- ^ / ^/dt dt' ^\{t)0^(t,£)^4^£\ (A.9) 

where we have made use of the reality of <I>^(x, t) to set = and defined 

= I \ . : : . .^ i Kt-^\ (a.io) 



2iB^b{t - to) -[h+^k- 2^fc^(* - *o) 

Afc, being the Fourier modes of the kernels above and the (potentially time-dependent) 
frequency. 

The Green's function Qk{t,t') is then defined via 

J dr Ok{t,T)gk{T,t') = -i6{t-t')h, (A.ll) 
where I2 = diag(l, 1) and we write 

G,(tt')-(^^''^'^''^^^"^''''A (A12) 

The off-diagonal terms in the kinetic operator in eq.(A.lO) offer a potential obstruction to 
the construction of the Green's function Qk{T,t'). We will argue, however, that we can make 
a consistent change of basis in which the kinetic operator is diagonal. We start by writing 
Ok{t,t') = 0^^'^\t,t') + Mk{t;to)5{t - t') where 

0(1,1).^. ( (^+-^2A.^5(t-to)) 
' ' \ -[^+ul-2AkR5{t-to) 

V -Dfc ^kl 
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where A^rj are the real and imaginary parts of the kernel respectively. Doing this 
extracts from the kinetic operator a piece proportional to rj = diag(l,— 1) and the term 
in the action involving this part of the kinetic operator is then invariant under a time- 
independent SU{1,1) rotation of ^^(t)- ^^(0 = ^k^^{t)- The full kinetic operator trans- 
forms as ©^(i, t') = {U^^y 0^{t,t')U^^; using this together with eq.(A.ll) relates the Green's 
functions as Gk{t,t') = Uu Qk{t,t') ul If we can choose the SU{1, 1) transformation Uk so as 
to diagonalize A^fe(t;io)) then we can focus on finding the Green's function in the diagonal 
case and transform back to the unrotated one in the end. 

To show that such an SU (1,1) transformation exists, consider taking Uk to be of the 
form 

^ / cosh0 — sinh^i (A 15) 

I — sinh 9 cosh 9 j 

This is a real SU{1, 1) transformation, and if we choose 

cosh 26 = , , sinh 26 = ^ , (A.16) 

^^h - Bl ^Al - Bl 

we diagonalize 7W/fc(t;to) to be proportional to the identity, 

{U^^y Mk{t-M U^^ = -2i^Al,-Bl 6{t - to) h. (A.17) 

The kinetic operator now becomes 

€j+ujI- 2Ckd(t - to)] 
OUM')= r ' ^ o ..... . (A.18) 




A^+ul-2Cl5it-to): 



with 



Ck = AkR + i\/Alj-Bl (A.19) 



For future reference, there is another way to parameterize the initial state that relates the 
kernels directly to various correlation functions. Let's write the kernels as [42] 

-^Ak = '-^^-^f, Bk='-^. (A.20) 

This parameterization is useful since is the two-point function (<I>g<I>_^)(to), ^.k^k is the 
(symmetrized) correlator between the field and its conjugate momentum 7r_^ at the ini- 
tial time, while the combination -|- o"|/4^^ is the momentum- momentum correlator. The 
parameter Cfc is a measure of how mixed the state is: Tr[/9^(to)] = (1/ak) < 1- Using this 
parameterization we find that 

2 



Ah-Bl={^J >0, (A.21) 
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so that the square root appearing in Ck is manifestly real. The transformation matrix Uk 
becomes 



and 

Ck = f+^^■ (A.23) 

We see then that we can reduce the problem of finding a Green's function with an arbitrary 
Gaussian initial density matrix to one where the mixing terms are absent. We focus on this 
case in the next subsection. 

A.l Green's functions for Gaussian pure states 

The Green's functions equations for the case with no mixing are 

-i5{t-t'), (A.24) 

0, (A.25) 

0, (A.26) 

i6{t-t'). (A.27) 

From these equations we see that setting G^~{t,t') = G^~'"(t, t')*, G^^{t,t') = G~l~{t,t')* is 
a consistent ansatz, so we can focus on calculating G^^ {t,t'), G^~ {t,t'). We set 

Gt-'M = f>(t)f<{t')e{t-t') + f<{t)f>{t')e{t'-t), (A.28) 

Gt-M = f<it)f>{t'). (A.29) 

In order that these Green's functions actually satisfy eqs.(A.24-A.27), we have to demand 
that 



d2 
dt2 




2GkS{t 


-io)^ 


\GV{t,t') 


d2 
dt2 




2Gk5{t 


-io)^ 


\Gt'M 


d2 
dt2 




2Gl5{t 


-io)^ 


\G,Ht,t') 


d2 
dt2 




2Gl5{t 


-to)) 


\Gr{t,t') 



+ u;l- 2Gk5{t - to) ] f<{t) = 0, (A.30) 



j2 
di2 

W[f^J^] =i, (A.31) 

f^imit) = f^*it)f^*it), (A.32) 



where W [fk:fk] Wronskian between and and the last two conditions are 

imposed only for t > to; the last condition enforces the requirement that the equal-time 
Green's functions all be equal to each other and real. 
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Let {hk{t)} be a complete set of mode functions satisfying 

^ + ui^hk{t) = 0, (A.33) 

hk{to) = hko, (A.34) 
W[hk,hl] = i. (A.35) 

The modes of interest can then be expanded as 

= a|/ifc(t) + ^^hUt) - 2f}{to)Ckgk{t)e{to - t), (A.36) 

where the last term will be used to match to the delta function term in eq.(A.30). 

The modes in eq.(A.36) must solve the full differential equation in eq.(A.30). In partic- 
ular, we have to match the discontinuity coming from the delta function term which leads 
to 



^ + ^l) 9k{t) = 0, gk{to) = 0, ^9k{t) 



= 1, (A.37) 

t=to 



where the initial conditions on gk{t) come from ensuring the absence of terms proportional 
to derivatives of 5{t — to) and canceling off the term proportional to 6{t — to) in eq.(A.30). 
This initial value problem can be solved for in terms of the modes {hk{t)} as 

gk{t) = i{hk{t)hl^-hl{t)hko), (A.38) 

where we have made use of the fact that the Wronskian W [hk, h^] is constant in time. 

The reality condition fk{t)fkii) = fk*i^)fk*i^) ^iH be satisfied if we set f>{t) = 
f^*{t) for t > to, which in turn implies = /3^*, = a^*- The Wronskian condition 
W [f^, f^] = i together with that for the modes {/ifc(t)} then implies 

W = {a>/3<-P>a<)W[hk,hl] \a>f - |/3>|' = 1. (A.39) 

This shows that for t > to, the modes f^{t) are Bogoliubov transforms of {/ifc(t), as 
might have been expected. 

Finally, we can solve for the coefficients a^, in terms of f^{to), /^(^o) ^ind the initial 
values of the modes {hk{t)} and their derivatives, 

«l = i (hlJ^ito) - (hl{to) - CkKo) ffito)) , (A.40) 

/3| = i [-hko§{to) + (Afc(to) - Ckho) /|(io)) • (A.41) 

It is interesting to note that while continuity of the modes enforces f^{to) = fk*{to)i this 
is not true of the time derivatives of these modes. In fact we can see that demanding that 
= (3^* leads to a statement about the discontinuity in these derivatives, 

h (*oj Ik [W 
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The last step in our calculation is to undo the SU{1, 1) transformation that brought 
us from the mixed state to the pure one and see how the Green's function transforms. As 
described in the previous subsection, the relevant transformation is 

gkit,t') = u^^ g'j,it,t') , (A.43) 

giving us the Green's function Qi^[t,t') for a mixed state in terms of the pure state Green's 
function Q'^^(t,t'). Using the parameterization in eq.(A.20) as well as eq.(A.22) we find that 

Gfit^t') = G',^\t,t') + {G',^{t,t') + G'+-(t,t')) (A.44) 

for a,b = ziz. For equal times we have 

Gt{t,t) = akf>{t)f<{t) = a, \f>{t)f . (A.45) 
Thus the equal-time functions G'^{t,t) in the mixed basis are all equal and real, as expected. 
A. 2 Application to the EFT of inflation 

The entire formalism described in this Appendix can be applied to the case of the vr field 
describing inflaton fluctuations in the decoupling limit. To do this, we perform the usual 
change to conformal time r] and conformally rescale vr so as to generate a kinetic operator of 
the form in eq.(A.lO). Thus, define the field % via 

X{x,ri) = V2M'^a{ri)TT{x,r]), (A.46) 

with = eAfp[ff^/Cg. The kernels in the initial state density matrix can now be defined in 
the X basis. In the de Sitter limit, we can take the modes to be the usual Hankel functions 
corresponding to the Bunch-Davies vacuum. 



2 3/2V s u \^ ^^^^ 

We are eventually interested in the Green's function for the gauge-invariant curvature 
perturbation C, = —Hn. This is given by 

where we have used eq.(A.46) to transform between the x a-iid vr fields. 
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